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ABSTRACT 

A Reynolds stress turbulence model has been implemented in the COMMIX code, together with 
transport equations describing turbulent heat fluxes, variance of temperature fluctuations, and dissipation of 
turbulence kinetic energy. The model has been verified partially by simulating homogeneous turbulent 
shear flow, and stable and unstable stratified shear flows with strong buoyancy-suppressing or enhancing 
turbulence. This article outlines the model, explains the verifications performed thus far, and discusses 
potential applications of the COMMIX-RSM code in several domains, including, but not limited to, 
analysis of thermal striping in engineering systems, simulation of turbulence in combustors, and predictions 
of bubbly and particulate flows. 


INTRODUCTION 

The computer code COMMDC-1C [ 1 ] describes single-phase, three-dimensional transient thermo- 
fluid dynamic problems in engineering systems. This code has provided the framework for the extension of 
the standard k-e turbulence model to a more sophisticated model based on 1 1 transport equations for the 
dependent variables describing turbulence: six transport equations for the components of the Reynolds 
stress tensor, three equations for scalar turbulent heat fluxes, one equation for variance of temperature 
fluctuations, and one equation for dissipation of turbulence kinetic energy. In isothermal calculations, the 
system reduces to seven transport equations. The complete Reynolds stress model (RSM) is presented in 
detail in the next section. A numerical verification of the reduced isothermal system has been made for 
homogeneous turbulence sustained by a uniform shear in the mean flow. In that case, the governing 
equations reduce to a system of ordinary differential equations that can be integrated directly, thus providing 
numerical verification of the code. Moreover, experimental data are available for comparison. The 
subsequent sections explain ongoing applications of the RSM to stratified flows, which also provide further 
code verification, and potential applications to several domains of interest for industrial applications where 
turbulence is highly anisotropic and the standard k-e model therefore performs poorly. The envisaged 
domains of application include analysis of thermal striping, simulation of turbulence in combustors, and in 
multicomponent flows. 

GOVERNING EQUATIONS FOR REYNOLDS STRESS MODEL 


The transport equations for scalar heat fluxes are 
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The transport equation for variance of temperature fluctuations is 
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The transport equations for Reynolds stresses are 
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The turbulence kinetic energy is k = u t u, / 2 . The transport equation for the dissipation of turbulence 
kinetic energy e is 
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APPLICATIONS OF REYNOLDS STRESS MODEL 

In this section, we summarize ongoing and potential applications of the RSM in the COMMIX code. 
The applications involve, but are not limited to, the following fields: 

• Stable and unstable stratifications of fluids in engineering systems. 

• Thermal striping. 

• Turbulence in combustors and in chemically reacting flows. 

• Bubbly and particulate flows. 

Stable and Unstable Stratifications o f Fluids in Engineering Systems 

In several engineering systems, including advanced nuclear reactors of the pool type, it is important to 
analyze the effect of fluid stratifications on turbulence and to predict under which operating conditions 
natural circulation flows can become established. Because turbulence in stratified flows is highly 
anisotropic, the standard k-e model (which assumes isotropy) performs poorly and must be replaced by 
more sophisticated models, such as the RSM. This domain of applications is exemplified by the two 
following test cases, which provide further verification of the RSM. 

Stratified Shear Flow 


An experiment designed to study the effect of buoyancy on turbulent mixing is the case of a horizontal 
shear flow [5] as sketched in Fig. 1. This figure shows qualitatively the spreading of the region where two 
fluid streams mix upon entering the domain with different temperatures and velocities. The influence of 
buoyancy is measured by the reduced Froude number 
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where the subscript z refers to the vertical component of the gravity acceleration. In stable stratification, 
when the hot fluid enters at the top, gravity forces oppose the diffusive character of turbulence, which tends 
to mix the fluids despite their density differences. The smaller the Froude number, the more stable the 
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Fig. 1. Sketch of stratified shear flow 
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The numerical values of the coefficients used in the above equations are as follows: c s + = 0.07, c^ = 
3.1, c 2 * = 0.4, c 3 * = 0.5, ci* = 0.5, c* = 0.13, R = 0.5, c k = 0.09, C1 = 2.8, c 2 = 0.47, c 3 = 0.4, c le = 
1.44, C2e = 1.92, c 3e = 0.8, c e = 0.15. Further details about the turbulence model and the boundary 
conditions used for all transport equations are given in Refs. [2] and [3]. 

VERIFICATION OF REYNOLDS STRESS MODEL 


The RSM has been verified for isothermal homogeneous turbulence maintained by a uniform shear in 
the mean flow. Homogeneity requires an infinite spatial field that can be simulated numerically by a finite 
domain with slip-free boundary conditions. Experiments on nearly homogeneous turbulent shear flow 
have been made in a wind tunnel of height h = 30.48 cm and length L = 10.5 h [4]. Turbulence generated 
by grids at the tunnel entrance becomes almost homogeneous in the center of the tunnel at some distance 
from the inlet. Measurements were taken in the region defined by 8.5 h < xj ^ 10.5 h, xj being a 
coordinate along the horizontal wind tunnel axis. Let x 2 and x 3 be coordinates along a vertical and 
transverse axis of the tunnel, respectively. In homogeneous turbulence, the mean velocity components of 
the flow satisfy the conditions U 1 =U 1 (x 2 );U 2 =U 3 =0; 8U 1 /3x 2 = constant (= 1 3 in the test case we 
refer to); U 2 (h/2) = U c = 12.4 m/s (center velocity). Under these assumptions, the seven transport 
equations describing the isothermal flow (Eqs. 6 and 7), reduce to a system of ordinary differential 
equations that can be integrated directly with the Runge— Kutta algorithm. We made numerical calculations 
of homogeneous turbulence with the RSM in the COMMIX code and with the Runge-Kutta method in an 
independent program; we then compared the results with experimental values at xi/h = 10 reported in 
Ref. [4]. A comparison of computed and experimental results is given in Table I. Under the assumptions 
made, uju^ and u^ are negligibly small. The computed results are much closer to each other than to the 
experimental values, due to the difficulty of realizing in practice the theoretical conditions of homogeneous 
turbulence. Theoretically, it must be u 2 u 2 = u 3 u 3 , a condition satisfied by the results of the calculations 
but not strictly verified in the experiment. This comparison provides satisfactory verification of the RSM 
under the given conditions. 

Table I. Comparison of experimental and computed turbulence values in homogeneous shear flow 

Turbulence Value Ref [4] Runge-Kutta COMMIX— RSM 

Vk/U^ 0.0185 0.0139 0.0147 

t/u^/u,. 0.0178 0.0140 0.0161 

V^/U c 0.0132 0.0097 0.0093 

0.0141 0.0097 0.0093 

■J|u 1 u 2 |/u c 0.0104 0.0101 0.0100 
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stratification. Below some threshold of the Froude number, turbulence is completely inhibited by the 
density gradients. In unstable stratification, when the hot fluid enters at the bottom, the effects of gravity 
and turbulence combine, forcing a region of highly effective mixing between the hot and cold fluids. 

In this study, we used 100 meshes (Ax = 0.03 m) in the axial direction along the channel length and 
30 meshes (Az = 0.006667 m) in the transverse direction. The maximum length (x max ) in the x direction is 
3 m and the distance (2h) between the two plates is 0.2 m. Therefore, the length-to-height ratio (x m ax/h) 
is 30. We present results of a test case with unstable stratification and Froude number 0.9, in which hot 
water is injected at the lower half of the test section, with Ui = U 2/2 = 0.223 m/s, Ti = 128.7°C, T 2 = 
20°C. 


Figure 2 shows the temperature distributions at location x/h = 20. With strong mixing, the 
temperature distributions in hot and cold fluids reach almost the same values. Compared with the k-e 
turbulence model, the RSM can produce results closer to the experimental results by accounting for 
anisotropy of natural circulation and turbulence transport. 

Thermal Stratification in a Rectangular Cavity 

Thermal stratification experiments have been made in a rectangular cavity filled with liquid sodium 
and simulating the hot plenum of a liquid-metal fast breeder reactor under steady-state and transient 
conditions [6]. The cavity communicates at the bottom with a rectangular channel (Fig. 3). A forced flow 
through the channel induces recirculation in the cavity. In steady-state experiments, temperature differences 
were induced by heating one wall of the cavity; in transient tests, by decreasing the temperature of the inlet 
flow. 


We simulated one steady-state experiment (referred to as PI in Ref. [6]). In this test, the inlet flow 
temperature was 300°C, while the right-hand wall of the cavity was kept at 316°C. Measured and computed 
vertical normalized temperature profiles in the cavity 200 mm from the left wall are given in Fig. 4. As in 
the previous test cases, the k-e model performs poorly, while the RSM results agree fairly well with the 
experimental ones. 

Thermal Striping 

Thermal striping is characterized by random temperature fluctuations in regions where flows with 
different temperatures mix. The temperature fluctuations, with frequency components of 2 to 20 Hz, induce 
cycle fatigue in structural materials. The problem of analyzing thermal striping and minimizing its impact 
upon structures is of concern in both nuclear and conventional industries. 

Combined buoyancy and striping phenomena have been investigated at ANL in piping, plena, heat 
exchangers, and steam generators. The distribution of temperature and velocity fields at the point where a 
horizontal pipe conveying water enters a plenum has been investigated as a function of the temperature 
difference between pipe inlet and bulk plenum temperature. In particular, the experimental program aimed 
at investigating (a) the length and height of the stratified recirculating flow carrying fluid from the plenum 
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Temperature Distribution (°C) 


Fig. 2. Temperature distributions for unstable stratified shear flow at x/h = 20 



Normalized Mean Temperature 

Fig. 3. Sketch of thermal cavity Fig. 4. Measured and computed normalized 

temperature distributions in the 
cavity, 200 mm from the left wail 

into the pipe; (b) the temperature fluctuations at the pipe walls, hence the potential damage due to thermal 
striping; and (c) the interaction between the thermal plume and the fluid plenum, hence the buoyancy- 
generated large-scale eddies in the plenum and the time evolution of plenum temperatures. 
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Numerical simulations of the flow and temperature fields at the junction between a horizontal pipe 
conveying hot fluid and a plenum initially filled with cold fluid have been made with the k-e turbulence 
model [7]. These computations failed to reproduce correctly the dynamics of the thermal plume in the 
plenum and the length of cold fluid penetration into the pipe. The latter was underestimated by a factor of 3. 
Thus, thermal striping analysis in the pipe could only be inferred. These discrepancies are due to the 
approximate representation of buoyancy forces and to the assumption of isotropy of turbulence inherent in 
the k-e model. Both of these limitations are eliminated in the RSM, which is more appropriate for detailed 
analyzis of this and related problems. Numerical simulation of this test case with the RSM is undertway at 
ANL. 

Turbulence in Combustors and in Chemically Reacting Flows 

Advanced combustion technologies are being developed with the goal of achieving higher overall 
system efficiencies and reduced environmental loading of air and water with solid pollutants. Further 
developments tend toward systems combining a pressurized, high-intensity combustor (which produces 
higher entrainment of the ambient fluid and fast mixing) with gas or steam turbines to optimize cycle 
efficiency and minimize pollutant formation. 

Computation of swirling turbulent flows in which turbulence is nonhomogeneous and anisotropic is a 
problem of major concern. It has been suggested that the assumptions leading to the formulation of the k-e 
model are inadequate for the simulation of turbulence in highly swirling flows [8]. Traditionally, designers 
have relied on experiments to derive empirical correlations describing heat and mass transfer. The validity 
of these correlations is generally limited to the experimental conditions for which they have been obtained. 
This traditional approach can now be supplemented with computational turbulence modeling (e.g., 
COMMIX-RSM), which gives the designer the predictive capability to change parameters and choose the 
conditions that lead to improved combustors with higher efficiency and less pollution. 

From the standpoint of numerical simulation of combustion involving several chemical components, 
the state-of-the-art code is KIVA-II [9] ; this is the last of a series of programs developed for, but not 
limited to, applications to internal combustion engines. This code is applicable to a variety of multi- 
dimensional fluid-dynamic problems with or without chemical reactions. Turbulent flows are computed 
with either a standard version of the k-e model or a subgrid scale (SGS) turbulence model in which large- 
scale eddies are resolved while small scale turbulence is modeled. These models, however, do not account 
for the presence of particles (e.g., coal particles in coal-fired diesel engines) in the carrying fluid and for the 
damping of the turbulence kinetic energy due to the presence of particles. Reynolds stress models can be 
applied to multifield turbulent transport by extending the procedure used for one-fluid turbulence to the 
several fields and taking into account their interactions [10]. The potential improvement of multiphase 
turbulence modeling provided by multifield RSM is worthy of consideration, at high computational costs, 
however. 

Bubbly and Particulate Flows 

Numerical simulations of bubbly and particulate flows in pipes are of interest for a variety of 
applications ranging from analysis of blood vessels to material transport by slurry flows. Numerical 
analysis aims at predicting pressure drops, relative velocity between continuous and particulate phases, and 
particle concentrations across flow stream lines. Turbulence predictions in the continuous phase have been 
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made [ 1 1] with a standard k-€ turbulence model and additional terms that account for the interfacial 
turbulent momentum transfer. Numerical analysis of turbulent bubbly flow and particulate flow is being 
carried out at ANL with the same approach as in Ref. [11]. Preliminary results for bubbly flows have been 
documented in Ref. [3], The assumption of isotropic turbulence inherent in the k-e turbulence model is 
however not satisfied in these cases. An extension of the RSM to two phases has been proposed in Ref. 
[12] and is being planned at ANL, in a way suitable for the numerical simulation of multicomponent flows. 

CONCLUSIONS 

In all cases where turbulence is anisotropic, and especially when enhancement or suppression of 
turbulence intensity due to buoyancy forces plays a dominant role, the RSM allows more realistic numerical 
predictions than the standard k-e model, which is based on the inherent assumption of isotropic turbulence. 
The RSM can be applied to a variety of domains of interest for technological applications. Its extension to 
multicomponent flows is possible, at the expense, however, of high computational costs. To cope with 
these, more advanced solution algorithms for the transport equations and code parallelization are envisaged. 
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NOMENCLATURE 


c p Specific heat (J/kg-K) 

F r Froude number 

Gk Production or suppression of turbulence kinetic energy due to buoyancy (J/s-m^) 
Gj* Buoyancy production in scalar flux equations (m-K/s 2 ) 

g ( = 2 )’ one_half of variance of temperature fluctuations (K 2 ) 

g Gravity acceleration (m/s 2 ) 
k Turbulence kinetic energy (m 2 /s 2 ) 

L Length scale (m) 

Pk Mean shear production in k and e equations (J/s-m^) 

Pi+ Mean field production in scalar flux equations (m-K/s 2 ) 

Rf (= - Gk/Pk). Richardson number 
T Temperature (K) 

t Time(s) 
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U Mean flow velocity (m/s) 
u Fluctuation of velocity (m/s) 
u i$ Scalar heat flux (m-K/s) 
u i u j Reynolds stress (m^/s^) 

xj Coordinate direction (m) 

Greek 

p «-( — ), volume expansion coefficient at constant pressure (K - * ) 

L pvaryj 

5ij Kronecker delta 

e Dissipation of turbulence kinetic energy (W/kg) 

X Thermal conductivity (W/m-K) 
p Dynamic viscosity (kg/m-s) 
v Kinematic viscosity (m^/s) 

71 ^ Pressure-scalar gradient correlation in scalar heat flux equations (m-K/s^) 
p Density (kg/m^) 

4» Temperature fluctuation (K) 

$ 2 3 4 V ariance of temperature fluctuations (K^) 

Indices 

l Laminar 
n Normal to wall 
t Turbulent 
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